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A key ingredient of cavity quantum- 
electrodynamics (QED) is the coupling between 
the discrete energy levels of an atom and photons 
in a single-mode cavity. The addition of periodic 
ultra-short laser pulses allows one to use such 
a system as a source of single photons; a vital 
ingredient in quantum information and optical 
computing schemes. Here, we analyze and 
"time-adjust" the photon-counting statistics of 
such a single-photon source, and show that the 
photon statistics can be described by a simple 
'transport-like' non-equilibrium model. We then 
show that there is a one-to-one correspondence 
of this model to that of non-equilibrium trans- 
port of electrons through a double quantum 
dot nanostructure. Then we prove that the 
statistics of the tunnelling electrons is equivalent 
to the statistics of the emitted photons. This 
represents a unification of the fields of pho- 
ton counting statistics and electron transport 
statistics. This correspondence empowers us 
to adapt several tools previously used for de- 
tecting quantum behavior in electron transport 
systems (e.g., super- Poissonian shot noise, and 
an extension of the Leggett-Garg inequality) to 
single-photon-source experiments. 

Cavity QED studies the interaction between a two- 
level atom and a single-mode cavity (see, e.g., refs. 1-4). 
Vacuum Rabi oscillations, the coherent excitation trans- 
fer between atoms and cavity photons, can occur if the 
atom-photon coupling strength g overwhelms both the 
loss rate (k) of the cavity photons and the emission rate 
(7) into other modes, as shown schematically in Fig. 1. 
To observe such a quantum oscillation, a velocity-selected 
atomic beam is passed through an open Fabry-Perot res- 
onator to control the interaction time ti. The probabil- 
ity P e {ti) that the atom remains in the excited state |e) 
at time ti is written 1-3 as P e = (1 + cos2gti)/2. This 
coherent coupling can be observed by examining the so- 
called vacuum Rabi splitting (VRS) in the transmission 
spectrum of the cavity. Clear evidence of VRS has been 
demonstrated not only in atomic systems 2-4 , but also 
in semiconductor self- assembled quantum dots 5 (QD) and 
circuit QED 6-8 systems. 

In several recent experiments (see, e.g., refs. 1-4), a 
cavity QED system was used as a source of single pho- 
tons by deterministically exciting the atom via periodic 



ultra-short laser pulses. Normally, one interprets the to- 
tal photon statistics from such an experiment as the en- 
semble average of a single event: the atom is excited at 
t = 0, interacts with the cavity, and eventually the cavity 
photon is emitted at some later time. All of the recorded 
single-photon detection events are then combined to give 
the ensemble average of this single situation. 

Here, we propose a simple alternative method of an- 
alyzing the photon detection events of this kind of 
experiments 1-4 that we term "time-adjusted photon 
counting" . We show that the photon emission spectrum 
can then be modelled via a Markovian master equation 
which has a one-to-one correspondence to a well stud- 
ied model of a double quantum dot (DQD) in the large 
bias, Coulomb-blockade regime. This allows us to reinter- 
pret data from existing (and future) cavity QED single- 
photon-source experiments as a continuous "transport- 
like" phenomenon, unifying photon and electron statis- 
tics. A summary of this correspondence can be found in 
Table I. 

Double quantum dots are artificial atoms in a 
solid 9-11 . A variety of powerful tools have been 
developed to study transport through such de- 
vices. These tools have revealed unique features like 
Coulomb-blockade 12 , the Kondo effect 12 , and coher- 
ent oscillations 13,14 . Our main result here is that the 
electron- counting statistics developed for the DQD model 
(e.g., current, current- noise, and higher-order cumulants) 
can be observed in existing photon- counting cavity QED 
experiments. We will use the analogy between these two 
apparently-unrelated systems to show that the photon 
statistics has a non-negative shot-noise feature, comple- 
mentary to their sub-Poissonian anti-bunching statistics, 
that indicates VRS. Moreover, we calculate the second- 
order correlation functions and show that these violate 
an extended form 15 of the Leggett-Garg inequality 16 . For 
completeness, we also consider violations of this inequal- 
ity by the non-adjusted statistics. 
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System 


Double Quantum Dot 


Cavity QED 


Carrier 


Electrons 


Photons 


Ground State 


\R) — electron in the right dot 


\g, 1} = ground state atom, 1 photon ) 


Excited State 


\L) = electron in the left dot 


|e, 0) = | excited atom, photons ) 


Energy difference AE 


El — Eb 


8/2 = {u-n)/2 


Rabi Rate 


Tunnelling amplitude T 


Atom-Photon coupling g 


Input Rate 


Tunnelling rate — > oo 


Laser pulses with time-adjusted shift 


Output Rate 


Tunnelling rate Tb 


Cavity loss rate k 


Quantum noise signature 


Super- Poissonian F e (Lu — >• 0) > 1 


Non- negative F p h(cj — > 0) > 


Extended LG inequality 


\2(I(t + r)I{t)) - (I(t + 2r)7(t))| < r«<7(t)) 


\2gW(t,t + T) -g {2) (t,t + 2r)\ < (a* (t)a(i)}" 1 



TABLE I: Comparison between the properties of the cavity QED system studied here and a double quantum dot. 



I. RESULTS 



B. Time-adjusted photon counting 



A. Standard photon-counting 

A way to produce 1 single photons from a cavity is the 
following: Ultra-short laser pulses with a given time con- 
stant are applied to the atom. The single-photon detector 
records the arrival time t n of a photon decaying out of 
the cavity with respect to each pulse, as shown Fig. 1. 
A normalized histogram of detection times reveals 1 pho- 
ton anti-bunching and oscillations due to the atom-cavity 
coupling. 

Neglecting the emission rate 7 into other modes, the 
normal cavity QED system can be described by the 
Markovian master equation (setting H = 1 throughout), 



p = W c [p] = -i[H c , p] + L c [p], 



where 



and 



L c [p] 



napa) — — [o)ap + pa) a] 



H c = va'a + —<j z + g((J-aJ + a + a). 



(1) 



(2) 



(3) 



Here to is the atomic level splitting, v is the cavity fre- 
quency, and k the cavity loss rate via which we acquire 
photons. Here, (a) denotes the creation (annihila- 
tion) operator of a cavity photon. The atomic opera- 
tors are defined as: a z = \e)(e\ — |<7)(<7|, cr_ = |ff)(e|, 
and (T + = |e)(g|, where |e) and \g) denote the excited 
and ground states, respectively. The atomic polarization 
decay 7 can be easily included in this analysis, but for 
simplicity we neglect it here. Furthermore, we omit vari- 
ations in the coupling strength g that can occur between 
each pulse. This could be an important factor, which can 
be dealt with by numerically integrating our final result 
over a Gaussian spread in g, or by including an additional 
dephasing term in the master equation. 



Now, rather than "collating" data in the above man- 
ner, as shown in Fig. 1(b), we propose to perform a 
time-adjusted analysis of the photon data, as shown in 
Fig. 1(c). Namely, the time between a photon-detection 
event and the next laser pulse (shown in green in the fig- 
ure) is eliminated, moving the time of each laser pulse 
to the time of the previous photon count; and any pe- 
riods of time with no-photon detection are eliminated 
from the data set. Thus, the system can then be viewed 
as one with instantaneous feedback that maintains one 
excitation in the combined atom-cavity basis. 

We show this rigorously by using the feedback for- 
malism in refs. 17-18 (refer to Methods for the detailed 
proof). Ultimately, this transforms the original equa- 
tion of motion [equations (1-3)] into an equation of mo- 
tion in the pseudo-spin two-state basis [defined as a z — 
|e,0)(e,0| — \g, l)(g, 1|, where e (g) represents the atom 
in its excited (ground) state, and 1 (0) represents a single 
(no) photon in the cavity]. The restricted master equa- 
tion (where feedback has been implicitly introduced, see 
Methods), is 



2 + 2 az + 9^x,P 



+ Kd+pa- - - 



<f- a+ p + p<f- <7+ 



(4) 



where S = uo — v is the detuning between the atom and 
the cavity. This restricted-basis equation of motion is 
equivalent to a two-level atom undergoing resonance flu- 
orescence in free space 19 . Here, the two-level atom is 
represented by the combined atom-cavity states \g, 1) and 
|e,0). The coherent input field is the natural atom-cavity 
interaction, and the time-adjusted photon counting gives 
an effective decay from \g, 1) to |e, 0), by eliminating 
the no-excitation state |<7,0). The time-adjusted data- 
set then represents a single trajectory in the ensemble 
described by this new equation of motion. 
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Laser pulse 

leakage 
photon 



Normal photon counting 




h h h 

Time-adjusted photon counting 



FIG. 1: (a) Schematic diagram of a QED system. Vac- 
uum Rabi oscillations can occur if the atom-photon coupling 
strength g overwhelms the loss rate « of cavity photons and 
the emission rate 7 into other modes, (b) Normal photon 
counting: Periodic ultra-short laser pulses excite the atom 
faster than all other time scales. The single-photon detector 
records the arrival time t n of a photon decaying out of the 
cavity with respect to each pulse. The grey area in (b) means 
that no photon is detected due to detector inefficiency, (c) 
Time-adjusted photon counting: The time (shown in green) 
between a photon detection event (in red) and the subsequent 
laser pulse is eliminated, moving the time of each laser pulse 
to the time of the previous photon count. Any periods (in 
grey) of no-photon detection are eliminated. In this manner 
the system has an effective instantaneous feedback: once a 
photon is detected, a laser pulse immediately drives the atom 
to the excited state, (d) The zero-frequency component of the 
power spectrum F(Q) is positive in the presence of coherent 
VRS (vacuum Rabi splitting), even if the coupling g is very 
small. In (d), brown is maximum, blue is minimum, and the 
black vertical lines show F(0) — 0. 



order correlation function 



(2) (q1(t)a1(t + T)a(t + T)a(t)) 

9 [ - t -- t + T >~ (at(t)a(i)) 2 W 



is found to be 



— [3k - 49 + 80e Qt - 4ae 2et ] , (6) 



89 



where a — 



^ + 9, 9 = ^k 2 /16-4 5 2 , and 5 = 
for convenience. It is important to point out that one 
cannot define a correct first-order correlation function 
G (1) (M + t) = (a f (t)a(t + t)) with this two-state model, 
unless one performs a full numerical simulation using a 
trace-preserving feedback operator (see Methods), or one 
retains the incoherent transition through the \g,0) state 
in the equation of motion. 



D. Analogy with electron transport 

As discussed earlier, our goal is to show that this simple 
model is equivalent to the electron transport through a 
solid-state double quantum dot, and then take advantage 
of common tools from transport theory. In the trans- 
port regime, a double quantum dot is connected to elec- 
tronic reservoirs with tunnelling rates Tl and Tr. If one 
assumes the strong Coulomb-blockade regime, i.e., the 
charging energy is much larger than other parameters, 
then one only needs to consider a single level in each dot. 
One can then define the three-state basis: \L), \R) and, 
|0) representing an electron in the left-dot, the right-dot, 
and the empty state, respectively. The DQD Hamilto- 
nian is written as 

H d = E L \L) (L\ + E R \R) (R\+T \L) (R\+T\R) (L\ , (7) 

where El (Eh) is the energy for the left (right)-dot 
level and T is the coherent tunnelling amplitude between 
them. 

The density matrix p(t) of the DQD satisfies 



dt 



p(t) = W d [ P (t)} = -i[Hd, P(t)} + L d [p(t)}. (8) 



The L d term in equation (8) contains the transport prop- 
erties and dissipation within the device, 



LMt)} 



[s L s L p(t) - 2s f L p{t)s L + p(t)s L s L '_ 



Second-order correlation function 



s i? s flP(*) ~ 2sitp(t)s]j + p(t)s' R s R ]l9) 



From this simple two-state model, the ensemble- 
averaged measurements of the photon output from the 
cavity can be easily calculated. In particular, the second- 



where s L - |0) (L\ , s[ = \L) (0| , s R = |0) (R\ , and s^ R = 
\R) (0|. One can calculate the current of electrons leaving 
the device using a current super-operator (e.g., for the 
junction on the right, and setting the electric charge e = 1 
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throughout) I R p(t) =T R \0) (R\ p(t) \R) (0|. The steady- 
state current of the right junction is then 20 ' 21 



(L) = Tr(W) 



r^T 2 



ry4 + e 2 + T 2 (2 + r K /r i )' 



(10) 



where e = E L — E R is the energy difference between the 
two dots. The shot noise S e (uj) of this device 22 ' 23 can 
also be easily obtained, 

/oo 
dre*- T <[5Jfl(t),W fl (* + T)] + )t_oo 
-OO 

+ 2</ a )5(r), (11) 

where the fluctuating right-junction current is SI R (t) = 
lB.{t) — (I s ), and the self-correlation term 2(I s )S(t) rep- 
resents the correlation of a tunnelling event with itself 20 . 
The shot noise (zero- frequency noise) is found to be 22 ' 23 , 



F e = S e (0)/2e(I s ) = 1 



8T 2 r L 



(12) 



Ae 2 {T R - r L ) + 3,T L T% 



r 3 + 8t r t 2 



-4r L e 2 + 4T 2 (r K + 2r L )] 2 



To reduce the double quantum dot problem, from a 
three-state basis to a two-state one, we take the limit 
of Tl S> T R ,T,e. This allows us to eliminate the |0) 
empty state, so that as an electron tunnels out of the 
right junction, an electron immediately tunnels through 
the left one; mimicking the effective two-level behavior 
of the cavity QED system. Interestingly, this limit gives 
exact results for all the stationary currents, but only re- 
turns the correct time-dependence for the right-junction 
current, in analogy with the inability for the restricted 
basis cavity-QED model to correctly construct In 
this case, the super-operator for the right-junction cur- 
rent becomes, I R p(t) = T R \L) (R\ p(t) \R) (L\. Similarly, 
the zero state is eliminated from the Ld term so that only 
one tunnelling rate, T^, remains. In this limited basis, 
equations (8) and (9) are equivalent to equation (4). Ta- 
ble I lists how the various parameters correspond to one 
another in the two different systems. 

Of particular interest to us is how the right-junction 
second-order current-correlation function in the large V l 
limit is equivalent to the second-order photon correlation 
function g( 2 > (t 7 t + r) discussed above. This is because (in 
the limit T L > T R ,T,e), 

(I R (t + r)I R (t)) = Tr[I R e w ^I R p(t)} 
= rlTT[\L)(R\{e w r\L)(R\p(t)\R)(L\}\R)(L\], (13) 

where each super-operator acts on those to the right. The 
corresponding correlation function for g^ 2 \t, t + r) for the 
cavity QED system, in the reduced basis we discussed 
earlier, is defined as 



(a f '(i)a f (t + r)a(t + r)a(i)) 
Tr[a_a+e w J^ + p(t)a_](14) 



where we have adopted the traditional input-output for- 
malism to define the photon intensity in terms of the 
internal pseudo-spin operators a. It is easy to see that 
the current-correlation measurement can be made equiv- 
alent to the second-order photon-intensity measurement 
simply by multiplying by a factor of k 2 . 

In summary, the ensemble-averaged photon statistics 
from a periodically-pulsed cavity QED system, following 
appropriate time-adjustments, has the same properties 
as the transport of electrons through a double quantum 
dot. As an example of the power of this apparently simple 
analogy, summarized in table I, we examine two tests for 
quantum behavior in double quantum dots (non-negative 
shot noise, and a special case of the Leggett-Garg inequal- 
ity) , and show how these two tests can be applied to the 
time-adjusted cavity QED system. 



E. Tests of quantum-ness 

1. Super-Poissonian shot noise 

It has been argued (and shown experimentally) that 
super-Poissonian shot noise can be observed in double 
quantum dots only if there is coherent quantum tun- 
nelling between the two dots 24,25 . If Fl is larger than 
r^, the second term in equation (12) becomes negative 
and produces a super-Poissonian value: S(0)/2e(I s ) > 1. 
In our analogy, Tl is always much larger than T R (it is 
effectively infinite). Using the correspondence between 
electron current and photon intensity, we can define an 
effective fluctuating photon-intensity noise-spectrum 



S* ph (w) = 2Re 



dre v 



where the effective photon current is (I p h) = K(a^ (t)a(t)) , 
and the subscript "ph" represents the photon analogy to 
typical electron transport measurements. In the pho- 
ton case there is no self-correlation term. We can easily 
calculate a photon current Fano factor using the same 
technique used for DQDs, and find 



Sph(O) 
2 (/ P h) 



8g 2 (3 K 2 -^ 2 ) 
(8.g 2 + k 2 + 5 2 ) 2 ' 



(15) 



From this equation one can easily see that super- 
Poissonian noise in the transport case corresponds to pos- 
itive noise in the photon case (F p h > 0). This can occur 
if 3k 2 < <5 2 , otherwise the 'shot noise' for photons is neg- 
ative. As mentioned above, in the DQD electron trans- 
port case, the super-Poissonian noise is only obtained 
for coherent coupling between the two dots. For clas- 
sical sequential tunnelling between two dots, the result 
turns out to be solely sub-Poissonian 24 ' 25 . For the time- 
adjusted single-photon cavity-QED system we consider 
here, coherent Rabi oscillations between the atom and 
cavity photon states produce a positive zero-frequency 
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component in the shot-noise spectrum. This is clearly 
indicated in Fig. 1(d), using parameters akin to those in 
the experiment in ref. 1. Even for g < k, the spectrum 
remains positive. Only for n 3> g does the whole range 
of Fph(O), as a function of S, become negative. 

In quantum optics, it is well known that a single- 
photon source has sub-Poissonian behavior in the second- 
order correlation function as the time r between mea- 
surements goes to zero: g^ 2 \t,t + t) t ^ — > 0. This is 
termed anti-bunching 26,27 . One can then call the light 
observed in this way as 'non-classical': it cannot be ob- 
tained from a thermal or coherent state source. Equation 
(13) implies a secondary criterium for the ' quantum- ness' 
of the observed light from such a single-photon source: if 
the correlation function is conditioned by quantum co- 
herent oscillations (vacuum Rabi splitting), one should 
see a positive value for the zero-frequency limit of the 
photon noise spectrum. 



k/2k = 2.7 MHz, 5 = 



y/2n = 3.0 MHz, k/2h = 2.7 MHz, g/2n 
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FIG. 2: A violation of the extended Leggett-Garg inequal- 
ity (15) for typical parameters in single-photon cavity QED 
experiments 2 ' 3 , using time-adjusted photon statistics. The 
parameters used here are: k/2tt =2.7 MHz, g/2n =10 MHz, 
and the variation of the detuning 5/2-k can be up to 20 MHz. 
The violation of the inequality is indicated by the grey "is- 
land" inside the black contour line. 



2. Extension of the Leggett-Garg Inequality 

To further clarify the quantum signatures in these 
correlation functions, we turn to an extension 15 of the 
Leggett-Garg inequality 16 . Recently, we derived our ex- 
tended inequality based on measurements of the fluctu- 
ating current through a DQD device 15 . Since the current 
is essentially an invasive measurement, we showed that 
the inequality only applies for systems with a state space 
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FIG. 3: (a) The non-adjusted photon-counting statistics also 
exhibit a violation of the Leggett-Garg inequality (15). The 
combined effect of the atom polarization decay and the cavity 
decay rates ( r y/2-K = 3 MHz and k/2-k = 2.7 MHz), prevent 
a violation being seen in the data in ref. 1. However, a slight 
decrease in the polarization decay rate 7; or an increase in 
the coupling strength g, combined with the ultra-short laser 
pulses in ref. 1, should reveal a violation of the inequality (15). 
(b) Using the same model (without time-adjusted counting), 
one can estimate the ratio between the coupling strength g 
and the dissipation (k, 7) needed to observe a violation of 
(15). This is approximately given by the vacuum Rabi split- 
ting VRS = ig/(n + 7) > 4. In (b) we have chosen k = 7. 



of three or less states, and under the assumption of irre- 
versible transport into the right reservoir (see Methods) . 
For the cavity QED system, these two assumptions be- 
come that of being in the single-excitation manifold, and 
the irreversible loss of photons once they leave the cavity, 
respectively. Thus, the current inequality of ref. 15, 

|Lj(t)| = \2{I(t + r)I(t)) (I(t + 2r)I(t))\ < T R {I(t)), 
becomes a second-order photon-correlation inequality, 



|L 9 (i)| = \2gW(t, t + r)- g™(t, t + 2r)\< (a\t)a(t)) \ 

(17) 

Superficially, this inequality is equivalent to the Leggett- 
Garg inequality in the stationary limit 16 . However, as 
mentioned above, the observation of a photon is an in- 
vasive measurement in terms of the cavity- atom state, 
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and thus, strictly speaking, this inequality no longer can 
be discussed in terms of distinguishing theories obeying 
macroscopic realism from quantum mechanics: the orig- 
inal goal of the work of Leggett and Garg. Here this in- 
equality only distinguishes quantum dynamics from those 
given by a classical rate equation. See the methods sec- 
tion, and ref. 15, for more details. Like the non-negative 
shot-noise feature, this inequality reveals a more nuanced 
way of understanding whether certain photon statistics 
have quantum characteristics, beyond those indicated by 
anti-bunching alone. 

In Fig. 2, we show how the violation of this inequal- 
ity occurs for a typical cavity QED experiment, using 1-4 
k/2tt =2.7 MHz, ■y/2ir = 3 MHz, g =10 MHz. As seen in 
Fig. 2, the violations of the inequality are easily observ- 
able and appear in a wide range of detuning 5. 



erties through a double quantum dot. This represents a 
unification of the fields of photon-counting statistics and 
electron-transport statistics. We then adapted several re- 
cent results from transport theory to describe or test the 
quantum nature of the photon-statistics being emitted 
from the cavity. 

We emphasize that not only the current noise, but also 
the higher-order moments 28 ' 29 can be examined with this 
time-adjusted scenario. Moreover, we point out that the 
same features could be observed in a circuit QED system, 
where the artificial atom (qubit) is periodically excited 
by some external means, and photons detected with a 
microwave photon counter 30 . 
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3. Extended inequality with standard photon statistics 

One can also apply the extended 15 Leggett-Garg 
inequality to the photon statistics without time- 
adjustment. In this histogram of the photon 
counts as a function of time (after the initial excitation 
of the atom) is equivalent to the second-order correlation 
function of the atom-cavity system with one excitation 
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system parameters. To give a more general bound for 
parameters which can cause a violation, in Fig. 3b we 
show the magnitude of the violation versus the cavity 
and atomic losses. This gives a bound on the Vacuum 
Rabi Splitting parameters needed to observe a violation 
of 2p/[(7 + k)/2] > 4. Many realizations of cavity QED 
have parameters which exceed this (see, e.g., ref. 5). We 
therefore think that this inequality (15) is a useful ad- 
dition to the toolbox one can use to test for quantum 
behavior in optical systems (see refs. 26,27 for reviews of 
other common tests). 
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II. DISCUSSION 



We have shown how a simple adjustment of the output 
photon-detection statistics of a periodically-excited cav- 
ity QED system can be described by a non-equilibrium 
model, with an exact analogy to electron transport prop- 
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III. METHODS 

Effective feedback formalism. First, we describe 
the measurement of a single photon (which has leaked 
from the cavity and is incident on a photodetector) as 
K{di)\ — \fdta. The complementary operator to this 
one, which is applied when no photon is detected dur- 



ing the duration dt, is K,(d£)2 — %/dia — (iH c + ^a l a)dt. 
The feedback evolution which is then applied to the sys- 
tem, given the observation of a photon, is described by 
0[...\ — |e)(g|...|<7)(e|; i.e. the atom is instantaneously 
and incoherently projected into its excited state. In gen- 
eral, this is not a trace-preserving evolution, as it is not a 
Liouvillian evolution, and is thus different from the class 
of feedback mechanisms described in rcfs. 17,18. Alterna- 
tively, one could assume a Liouvillian evolution accord- 
ing to the coherent dynamics of a laser pulse causing 
7r-pulse/Rabi oscillations of the atom from its ground to 
excited state, via an operator e~ %z , Z = [a x ,...]. How- 
ever, in the limit when this transition is faster than all 
other dynamics, and we are in the single-excitation man- 
ifold, then it is equivalent to O [...]. 

The unnormalized density matrix following the detec- 
tion of a photon at time t, and evolution due to feedback, 
is 

pi{t + dt) = nOap(t)a^dt. 

Our time-adjustment scheme implies that the time delay 
is far smaller (effectively zero) than the cavity decay time, 
and one can assume the fully Markovian master equation, 



P 



-i[H c , p] + KOapo) 



a'ap + pa) a 



If the feedback is via the operator O [...], this is not 
a trace-preserving equation of motion. However, if we 
restrict ourselves to the single (lowest) excitation man- 
ifold, only the |e,0), |g,l), and \g,0) states are im- 
portant. Then the feedback term becomes Oapa) — 
|e, 0}(g, 1|/?|<?, l)(e, 0|, and this truncated basis is trace- 
preserving, as the state |e, 1) is decoupled from the 
equation of motion, giving Trpapa^} = Tr[ apa 1 \ = 
{<!■ 1 I' .'/• 1 • 

Now the action of the instantaneous feedback is obvi- 
ous, such that the state \g, 0) in the photon-decay terms 
is decoupled from the single-excitation manifold. We can 
now write our equation of motion purely in the pseudo- 
spin two-state basis , denning a z — | e , 0) ( e , 1 — \g,l)(g,l\, 



P = -i[v + 5a z + ga x ,p\ 



Kcr + p(j- 



(J-U+p + /9(7-(7+ 



where 5 = lu — v is the detuning between the atom and 
the cavity. 

To achieve the above, we note the following points. 
First, the delay between measurement and feedback is 
instantaneous. In our case the delay is zero, as dictated 



by our adjustment of time intervals in the data set. Sec- 
ond, the feedback action (i.e., laser pulse) incoherently 
and instantaneously projects the system into the |e, 0) 
state. This has recently been achieved by Bochmann et 
al. 1 using ultra-short laser pulses. Third, the photon de- 
tection is here assumed to be 100% efficient. This can be 
effectively achieved by simply omitting the time intervals 
where no photons are detected. 

Derivation of the extended Leggett-Garg inequal- 
ity. In ref. 15, we extended the Leggett-Garg inequality 
to work under the conditions of invasive measurement, 
but with additional restrictions. We summarize and re- 
formulate the proof of that inequality here, but now using 
the language of cavity-QED. 

For the cavity-QED case, we posit that any photon- 
intensity measurements not conditioned by quantum dy- 
namics obey 



\L g (t,t+r)\ = \2g^(t,t+T)-g^(t,t+2r)\ < (a^t)a(t)} " 

In the language of an effective photon current, we can 
write this as 

|£/(M + r)| = (19) 
\2(I ph {t + T)I ph {t)) - (Jph(t + 2r)7 ph (t))| < K(/ ph ), 

where k is the rate of photons leaking from the cavity, 
I p h(t) = I p h(t = 0) and (I p h{t)} is the average photon 
current of the initial state. Hereafter we omit the t vari- 
able. In the master equation approach, the current oper- 
ator translates into a "jump" super-operator, and equa- 
tion (|19[) thus represents an inequality concerning transi- 
tions in the system, and not static properties. Thus it is 
obviously suitable for application to single-photon mea- 
surements, which give us information about a change in 
the state of the cavity-QED system. As described in 
the text, the photon current super-operator acts as be- 
fore / p h[p] = ncr + pcr~ , such that the average current is 

(^ph) = Tr |/ p hp| and the correlation function of interest 

is obtained as (I p h(T~)I p h) — Tr |/ p he £T ^phPo|- For our 
time-adjusted case, the stationary distribution is chosen 
as the initial state. For the non-time-adjusted photon 
statistics, one chooses po — \g, l)(l,g|. 

In these terms, the inequality expression can be written 

as L/(r) = Tr |l ph (2e £r - e 2Cr ) i ph p Q } . If the cavity- 
QED system contains no coherent quantum dynamics, 



Jph is the 3x3 matrix with elements I p h 



a/3 



where the indices = \g,0),V = \g,l),A= |e,0). Thus 
using the Chapman-Kolgomorov equation, we have 



Li(r) 



k 2 P v (Q) 



^vo (2 — Qoo — ^vv) — ^VA^AO 



Where fl represents the matrix elements of the propa- 
gator. For a general Markov stochastic matrix, fi, the 
maximum of Lj is 

max{£/} = 2k 2 P v (0). 
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However, the rate equation form S1(t) = exp Ct furnishes 
us with a further constraint. Maximizing Lj(t) with re- 
spect to time, from Li = and fl = £fi, we find that 
the maximum of Lj occurs when f2oo + &>vv — 1 arL d 
ttvo = 1, giving 

max{£/} = k 2 P p (0) = n(I ph ). 

This result relies on the form of the jump operator and 



the absence of re-absorbtion of photons by the cavity, 
i.e. C-po = 0. These requirements mean that we must 
always be in the single-excitation regime, and that once 
a photon leaves the cavity and is measured, it cannot 
return. Fortunately, this is implicit in the definition of 
destructive photon measurement. 
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